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ELECTROSTATIC  EQUATIONS  FOR  LARGE  SCALE  PLASMA  SIMULATION  STUDIES 


I .  INTRODUCTION 

In  recent  years,  global  numerical  simulation  of  plasma  dynamics  has 
become  an  area  of  active  research  interest.  A.  large  portion  of  this  effort 
has  been  based  on  the  equations  of  magnetohydrodynamics,  MHD.  The 
equations  of  MHD  include  a  plasma  continuity  equation,  a  momentum  equation, 
an  energy  equation.  Maxwell's  equations,  and  Ohm's  law.  These  equations 
allow  a  simulation  of  the  temporal  evolution  of  a  fluid  plasma  and  the 
electromagnetic  field.  Perpendicular  to  the  magnetic  field  the  simulations 
maintain  a  balance  between  plasma  forces  and  field  forces,  while  parallel 
to  the  field  the  balance  depends  only  on  fluid  quantities.  The  balance 
between  plasma  and  field  forces  is  responsible,  in  part,  for  the 
reliability  and  accuracy  of  the  simulation  results.  In  certain  problems  of 
interest  the  plasma  has  a  relatively  low  temperature  and  the  magnetic 
forces  dominate  the  motion  perpendicular  to  the  field.  This  situation 
leads  to  difficulties  in  accurately  computing  the  plasma  forces  and  therefore 
the  dynamics  of  the  plasma.  In  such  situations  the  difficulty  can  be 
alleviated  by  restricting  the  MHD  equations  to  the  so-called  electrostatic 
approximation.  In  the  electrostatic  approximation  one  retains  the  hydro- 
dynamic  equations  parallel  to  the  magnetic  field  but  perpendicular  to  the 
magnetic  field  balances  plasma  forces  against  J  x  B,  where  B  Is  a  known 
invariant  magnetic  field  and  J  is  calculated  from  the  difference  between 
electron  and  ion  velocities. 

The  purpose  of  this  report  is  to  derive  a  set  of  electrostatic 
equations  which  can  be  used  to  develop  a  simulation  code  for  large  scale  'j 

low  pressure  plasma  systems.  This  is  done  in  a  general  coordinate  I 
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system.  The  organization  of  the  paper  is  as  follows.  In  the  Section  II 
the  meaning  of  the  electrostatic  approximation  is  discussed.  In  Section 
III,  the  physical  equations  are  derived.  In  Section  IV,  the  electrostatic 
equations  are  written  for  a  general  coordinate  system.  In  Section  V,  the 
application  of  these  equations  to  a  plasma  system  in  a  dipole  coordinate 
system  appropriate  for  the  earth's  ionosphere-magnetosphere  is  presented. 
Finally,  in  Section  VI  we  briefly  summarize  our  results. 

II.  THE  ELECTROSTATIC  APPROXIMATION 

The  electrostatic  approximation  consists  of  four  assumptions: 

1)  All  forces  are  small  compared  to  the  magnetic  force  so  that 
3B/3t  ■  0  and  V  *  B  *  0.  In  order  for  B  not  to  change  in  time  it 
is  sufficient  and  necessary  that  the  electric  field,  is  the 
gradient  of  a  scalar,  i.e.,  E  *  -VHf . 

2)  The  current  J  is  given  explicitly  by  the  difference  in  the  ion 
and  electron  velocities,  i.e.,  J  -  en^V^  -  VgJ  .  The  requirement 
that  V* J  *  0  provides  an  equation  for  the  electric  potential  V  . 

3)  The  electrostatic  potential  is  constant  along  magnetic  field 
lines,  i.e.,  the  field  lines  are  equipotentials. 

4)  In  order  to  compute  the  velocities  perpendicular  to  the  magnetic 
field  lines  one  uses  the  approximation  JJ  +  x  ,?  ■  0.  This 
assumption  is  not  necessary  but  greatly  simplifies  the  analysis. 

Within  the  context  of  these  assumptions,  i.e.,  the  electrostatic 
approximation,  we  derive  an  appropriate  set  of  equations  for  performing 
plasma  simulations  in  a  general  coordinate  system. 


III.  DERIVATION  OF  EQUATIONS 


The  momentum  equations  for  the  electrons  and  ions  are,  respectively, 


-  —  i  E  +  V  x  Bl  -  v  (.V  -  V  j  -  v  I  V  -VI  -  g 

m  -e  env~e  "-n  ei  ~e  ~i  *6 


(1) 


^  u  +  h  *  s>  -  '’i.U  i  -  V  -  vi.Ui  -  ~V  ■  •*! 


(2) 


where  B  is  the  magnetic  induction,  E  is  the  electric  field,  V  ,  V  ,  and  V 

are  the  electron,  ion,  and  neutral  velocities,  respectively,  v  and 

are  the  electron-neutral  and  electron-ion  collision  frequencies,  v  and 

m 

v  are  the  ion-neutral  and  ion-electron  collision  frequencies, 
respectively,  and  ga  is  the  total  force  (including  inertia)  acting  on 
species  a. 

Equations  (1)  and  (2)  can  be  written  in  the  form 
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Eliminating  the  croaa  term,  from  Eqa.  (3)  and  (4)  one  finds  that 
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one  can  solve  for  V  and  V 
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We  now  derive  an  expression  for  the  current  J  using  Eqs.  (7)  and 
(8 )  .  We  note  that 
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tn  order  to  get  an  equation  for  the  potential  one  has  to  first  Isolate 
the  inertial  terms.  Neglecting  electron  inertia  and  defining 
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one  finds  that 
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We  rewrite  Eq.  (11)  using  the  following  notation 
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so  that  one  finally  obtains 
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The  total  time  derivative  d/dt  also  contains  all  coriolis  and  centrifugal 


terms. 


Using  7* J  ■  0  and  integrating  along  the  fields  with 


E  =  —  7Y  ;  V  -  (E  x  B )/ | B | - 


one  can  derive  the  following  two  dimensional  potential  equation  for  Y 


J dl  7*[op7Y+  (p/B)Tpld(VY/B)/dtj 


+  j dl  V*[{(-0 VY  +  (p/B)T  d(VY/B)/dtJ  x  B } /B J 
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This  equation  will  be  used  in  the  subsequent  analysis.  The  first 
integral  represents  the  Pedersen  conductance  owing  to  collisions  and  the 
capacitance  associated  with  plasma  inertia.  The  next  integral  represents 
the  analogous  Hall  terms.  The  first  two  terms  on  the  right  hand  side  are 
associated  with  collisional  coupling  to  the  neutral  gas.  Thev  are,  in  many 
instances,  the  main  driving  terms. 

One  should  remark  also  that  the  first  integral  in  Eq .  (13)  must  to  be 

larger  than  the  second  Integral  for  the  potential  equation  to  be  diagonally 

dominant.  By  looking  at  the  formulae  for  a  and  <j.  one  finds  that  the  field 

p  h 

line  integral  of  the  product  of  the  plasma  density  and  the  Larroor  frequency 
must  be  larger  than  the  comparable  integral  of  the  collision  frequency. 


IV .  GEOMETRY 

We  now  derive  an  expression  for  the  potential  equation  in  a 
generalized  geometry.  An  orthogonal  coordinate  system  (x^,X2»X3)  is 
defined  by  the  geometrical  factors  (h^.l^.h-j)  such  that  the  length  element 
ds  is  given  by 


2  2  2  2  2  2  2 
ds  -  hj  dxj  +  h2  d*2  +  h*  dx^ 
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The  covariant  derivative  of  the  velocity  is  given  by 
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Since  the  electrostatic 
one  of  the  coordinates  can 
follows  then  that 


approximation  assumes  a  vacuum  magnetic  field 
be  the  field  lines.  Here  x^  is  chosen.  I 
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which  means 


hlh2B3  "  fixl»x2') 


(20] 


is  a  function  of  x,  and  x0  only. 


with  the  definition  of 


dS3  ■  hj  dx^ 

the  potential  equation  (Eq.  (13))  can  be  written  explicitly  as 
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It  is  the  goal  to  separate  geometrical  quantities  from  fluid 


quantities.  To  this  end  we  define  normalized  quantities  at  a  point  x^«  We 


define  X^,  B®,  and  g®  as  follows: 
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The  only  Integrals  which  are  not  cast  in  the  standard  form  are  the 
transport  terms.  We  define 
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using  equation  (4). 

Furthermore,  we  note  that 
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We  define  Che  quantities 
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The  potential  equation  for  ¥  can  then  be  written  in  the  final  compact 
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V.  APPLICATION  TO  DIPOLE  COORDINATES 

For  the  electrostatic  approximation  we  assume  an  undisturbed 
geomagnetic  field  which  is  taken  to  be  a  pure  dipole  field.  Therefore,  we 
use  the  appropriate  dipole  coordinate  system  (n,$,c)  where 
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(50) 


where  r  is  the  radius,  <j>  the  magnetic  longitude,  0  the  magnetic 
colatitude.  An  orthogonal  coordinate  system  is  defined  by  the  geometrical 
h  factors,  which  we  derive  next. 

We  note  that 
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which  gives  the  Jacobian  D 
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The  inverse  matrix  is 


The  coefficient  for  the  gravity  is 
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sin6 

(1+3  cos20)1/2 


,02 
800(-r  J 


cos  9 


M  7  , 2.1/2 

[1  +  3  cos  9j 


(61) 


For  high  altitude  raagnetospheric  problems  one  can  set  vgn  *  0.  Then 

2 

with  *  eB/m^  ■  Q  X  and  b  ■  (v^n/Xjn^e  ■  ”  v/(fljX  J  and  ujje  ■  Qj  + 

n.  v.  one  finds 
ie  in 


J).  U), 

T  -  1  ie  . 

T  m 

,  2 
fl.ii).  +  v. 
i  ie  in 

(62) 

P  r  » 

0)  +  V 

ie  in 

pe 

2  ,  2 
il).  +  V. 

ie  in 

n.  v. 

T  —  i  in  . 

T 

2 

n,  v. 
ie  in 

(63) 

h  2  2  ’ 

0).  +  V. 

ie  in 

ne 

2  ^  2 
“ie+  vin 

S  -  pv,  T  ; 

P  in  P 

sh- 

-  pvinTh 

(64) 

Furthermore  we 

note  that 

Si"  s  "  P  7pi 

f  - 

~e 

-  -  Vp 

P  e 

(65) 

which  leads  to 

F  -  T  a  -  - 
-p  p  z  p 

lTp  ,pi+  V«)! 

F,  " 
-h 

Ti,r;l’h  'v  V'”.)- 

(66) 

We  choose  the  geomagnetic  equator  (a  •  0)  as  a  reference  plane.  At 
the  equator  one  has 

h1  -  l  ;  h2  -  n  ;  h^  -  l  ;  g  -  n  ;  b  «  A/n  (67) 
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The  final  equation  can  now  be  written  as 


n  3n  nn  3r^  2  34  34^ 


Ln  3n  v  nn  3n  St'  ^2  3$  *•  4>4>  3$  ZtJ1 


l  ,  3cn4>  32Y  3Cn$  32y  , 

+  n  TiTTt  3n  3t-* 


7T  tnlD!  +  T!  +  Ff  “  Dl!  +  T!?  +  FlJ] 

n  3n  1  v  ^  ♦  <P  n  n  $ 


~  ‘S'  I4  l-Dn  +  Tn  +  Fn  +  D*  -  TA  "  F*J  * 

n  3)  1  r)  n  n  $  $  * 


where  all  symbols  have  been  previously  defined. 


VI .  SUMMARY 


We  have  presented  a  set  of  equations  for  large  scale  simulation 


studies  of  plasma  systems  in  the  electrostatic  approximation.  The 


equations  are  written  for  a  generalized  geometry  as  well  as  for  a  dipole 


geometry  which  is  suitable  for  the  earth's  magnetosphere. 
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JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  2081Q 

01  CY  ATTN  DOCUMENT  LIBRARIAN 
01 CY  ATTN  THOMAS  POTEMRA 
01 CY  ATTN  JOHN  DASSOULAS 

KAMAN  SCIENCES  CORP 
P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 
01 CY  ATTN  T.  MEAGHER 

KAMAN  TEMPO-CENTER  FOR  ADVANCED 
STUDIES 

816  STATE  STREET  (P.O  DRAWER  00) 

SANTA  BARBARA,  CA  93102 
01 CY  ATTN  D  AS  I  AC 
01 CY  ATTN  WARREN  S.  KNAPP 
01 CY  ATTN  WILLIAM  MCNAMARA 
01 CY  ATTN  B.  GAMBILL 

LINKABIT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 
01 CY  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  8  SPACE  CO.,  INC 
P.O.  BOX  504 
SUNNYVALE,  CA  94088 
01 CY  ATTN  DEPT  60-12 
01 CY  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  8  SPACE  CO.,  INC. 
3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

01 CY  ATTN  MARTIN  WALT  DEPT  52-12 
01 CY  ATTN  W.L.  IMHOF  DEPT  52-12 
01 CY  ATTN  RICHARD  G.  JOHNSON 
DEPT  52-12 

01  CY  ATTN  J.B.  CLADIS  DEPT  52-12 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32305 

01 C Y  ATTN  R.  HEFFNER 


M.I.T.  LINCOLN  LABORATORY 
P.O.  BOX  73 
LEXINGTON,  MA  02173 

01 CY  ATTN  DAVID  M.  TOWLE 
01 CY  ATTN  L.  LOUGHLIN 
01 CY  ATTN  D.  CLARK 


MCDONNEL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
01 CY  ATTN  N.  HARRIS 
0 1 C  Y  ATTN  J.  HOULE 
0 1 C Y  ATTN  GEORGE  MROZ 
01 CY  ATTN  W.  OLSON 
0 1  C  Y  ATTN  R.W.  HALPRIN 
0 1 C  Y  ATTN  TECHNICAL 

LIBRARY  SERVICES 

HI  SSI  ON  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  93101 
01  CY  ATTN  P.  FISCHER 
01  C  Y  ATTN  W.F.  CREVIER 
01 CY  ATTN  STEVEN  L.  GUTSCHE 
01 CY  ATTN  R.  BOGUSCH 
01 CY  ATTN  R.  HENDRICK 
01 CY  ATTN  RALPH  KILB 
01 CY  ATTN  DAVE  SOULE 
01 CY  ATTN  F.  FAJEN 
01 CY  ATTN  H.  SCHEIBE 
0 1 C  Y  ATTN  CONRAD  L.  LONGHIRE 
01 CY  ATTN  8.  WHITE 
01 CY  ATTN  R.  STAGAT 

HISSION  RESEARCH  CORP. 

1720  RANDOLPH  ROAD,  S.E. 
ALBUQUERQUE,  NEW  HEXICO  87106 
01 CY  R.  STELLINGWERF 
01  C  Y  H.  ALHE 
01 CY  L.  WRIGHT 


PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

01  CY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHOTOHETRICS,  INC. 

4  ARROW  DRIVE 
WOBURN,  HA  01801 

01  CY  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAHICS,  INC. 

P.O.  BOX  3027 
BELLEVUE,  WA  98009 

0 1 C  Y  ATTN  E.J.  FREMOUW 

PHYSICAL  DYNAHICS,  INC. 

P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOHSON 

R  ft  D  ASSOCIATES 

P.O.  BOX  9695 

HARINA  DEL  REY,  CA  90291 


01  CY 

ATTN 

FORREST  GILMORE 

01  CY 

ATTN 

WILLIAM  B.  WRIGHT 

01  CY 

ATTN 

WILLIAM  J.  KARZAS 

0 1  C  Y 

ATTN 

H.  ORY 

01  CY 

ATTN 

C.  MACDONALD 

01  CY 

ATTN 

R.  TURCO 

01  CY 

ATTN 

L.  DeRAND 

01  CY 

ATTN 

W.  TSAI 

MITRE  CORPORATION,  THE 
P.O.  BOX  208 
BEDFORD,  MA  01730 

0 1 C  Y  ATTN  JOHN  MORGANSTERN 
01 C  Y  ATTN  G.  HARDING 
01 CY  ATTN  C.E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 
01 C  Y  ATTN  W.  HALL 
01 CY  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
12340  SANTA  MONICA  BLVD. 

LOS  ANGELES,  CA  90025 

01 C  Y  ATTN  E.C.  FIELD,  JR. 


RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
01 CY  ATTN  CULLEN  CRAIN 
01 CY  ATTN  ED  BEDROZIAN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  HA  01776 

0 1  C Y  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
330  WEST  42nd  STREET 
NEW  YORK,  NY  10036 

01  CY  ATTN  VINCE  TRAPANI 


SCIENCE  APPLICATIONS,  INC. 
1150  PROSPECT  PLAZA 


LA  JOLLA 

,  CA 

92037 

01  CY 

ATTN 

LEWIS  M. 

LINSON 

01  CY 

ATTN 

DANIEL  A. 

HAMLIN 

01  CY 

ATTN 

E.  FRIEMAN 

01  CY 

ATTN 

E.A.  STPAKER 

01  CY 

ATTN 

CURTIS  A. 

SMITH 

SCIENCE 

APPLI 

CATIONS, 

INC 

1710  GOODRIDGE  DR. 
MCLEAN,  VA  22102 
01 CY  J.  COCKAYNE 
01 CY  E.  HYMAN 


SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 


MENLO  1 

PARK, 

CA  94025 

01  CY 

ATTN 

J.  CASPER 

01  CY 

ATTN 

DONALD  NEIL SON 

01  CY 

ATTN 

ALAN  BURNS 

01  CY 

ATTN 

G.  SMITH 

01  CY 

ATTN 

R.  TSUNODA 

01  CY 

ATTN 

DAVID  A.  JOHNSON 

01  CY 

ATTN 

WALTER  G.  CHE5NUT 

01  CY 

ATTN 

CHARLES  L.  RINO 

01  CY 

ATTN 

WALTER  JAYE 

01  CY 

ATTN 

J.  VICKREY 

01  CY 

ATTN 

RAY  L.  LE AD A3R AND 

01  CY 

ATTN 

G.  CARPENTER 

01  CY 

ATTN 

G.  PRICE 

0 1  C  Y 

ATTN 

R.  LIVINGSTON 

01  CY 

ATTN 

V.  GONZALES 

01  CY 

ATTN 

D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

0 1 C Y  ATTN  W.P.  BOQUIST 

TOYON  RESEARCH  CO. 

P.O.  Box  6890 
SANTA  BARBARA,  CA  93111 
01 CY  ATTN  JOHN  ISE,  JR. 

01 CY  ATTN  JOEL  GARBARINO 

TRW  DEFENSE  8  SPACE  SYS  GROUP 


ONE  SPACE  PARK 


REDONDO 

BEACH 

* 

CA  90273 

01  CY 

ATTN 

R. 

K.  PLEauCH 

01  CY 

ATTN 

S. 

ALTSCHULER 

0 1  C  Y 

ATTN 

D. 

DEE 

01  CY 

ATTN 

D/ 

STOCKWELL 

SNTP / 1 575 

VISIDYNE 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MASS  01803 
0 1 C  Y  ATTN  W.  REIDY 
01 C  Y  ATTN  J.  CARPENTER 
01CY  ATTN  C.  HUMPHREY 

UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA  15213 

01 CY  ATTN:  N.  ZABUSKY 
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